The phase transition in the localized ferromagnet EuO probed by fiSH 
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We report results of muon spin rotation measurements performed on the ferromagnetic semicon- 
ductor EuO, which is one of the best approximations to a localized ferromagnet. We argue that 
implanted muons are sensitive to the internal field primarily through a combination of hyperfine 
and Lorentz fields. The temperature dependences of the internal field and the relaxation rate have 
been measured and are compared with previous theoretical predictions. 



Europium oxide (EuO) crystallises in the rock-salt 
structure and is a ferromagnetic semiconductor with a 
Curie temperature (Tq) of 69 K It shows a colossal 
magnetoresistance effect in its Eu-rich form 0, Q and 
the associated metal-insulator transition has been linked 
to the formation of bound magnetic polarons It is 
the only magnetic binary oxide known to be thermody- 
namically stable in contact with silicon [5|, and this, to- 
gether with a nearly 100% spin polarization of mobile 
electrons for carrier concentrations below half filling of 
the conduction band 0, 0] , means that it is thought to 
be highly relevant for spintronic applications The 
magnetic moments of the Eu 2+ (4f 7 , 8 S 7 / 2 ) ions result 
from 4f charge density which is nearly completely local- 
ized inside the filled 5s 2 5p 6 shells and there is negligible 
overlap. This makes EuO an excellent approximation to 
a Heisenberg ferromagnet although there is evidence 
of some momentum dependence in the exchange inter- 
actions (lo| . The ferromagnetic interactions are due to 
an indirect exchange mediated by electrons in anion va- 
lence bands l_l| , specifically virtual excitations of oxygen 
valence band p electrons into empty Eu 2+ (5d) conduc- 
tion bands and exchange interaction of the d electron (p 
hole) with the localized 4f electrons. EuO is one of a 
family of isostructural europium chalcogenides EuX with 
X=0, S, Se and Te. As the size of the chalcogen increases 
from O to Te, the lattice parameter a increases across 
the series: 5.14 A in EuO, 5.97 A in EuS, 6.20 A in EuSe 
and 6.60 A in EuTe Only EuO and EuS are fcrro- 

magnets and both the nearest-neighbor and next-nearest- 
neighbour exchange constants are larger in EuO com- 
pared to EuS @, In EuO, the saturation magnetiza- 
tion M s is given by gJ^ B (4/a 3 ), yielding ^i M s = 2.40T 
(using gj = 7). 

Neutron studies of EuO are significantly hindered by 
the strong absorption of neutrons by Eu, though the use 
of a thin-slab geometry and use of 153 Eu (which has the 
smaller absorption cross section of the two naturally oc- 
curring isotopes) has allowed a detailed study to be per- 
formed [9J. Nevertheless, a technique such as muon-spin 
rotation (^iSR) avoids these difficulties entirely, and in 
contrast to NMR there is no electric field gradient or 



quadrupolar contribution to the observed muon response, 
simplifying the analysis. In this paper we present the re- 
sults of /iSR experiments to study the magnetic order and 
the fluctuations in samples of EuO and use the results to 
compare with theoretical predictions EMI. 

Our /xSR experiments were carried out using the GPS 
instrument at the Swiss Muon Source (S/iS), Paul Scher- 
rer Institute (PSI) in Switzerland. In our zxSR exper- 
iment, spin polarised positive muons (p + , momentum 
28 MeV/c) were implanted into small crystals of EuO. 
The muons stop quickly (in < 10 -9 s), without signif- 
icant loss of spin-polarization. The observed quantity 
is then the time evolution of the average muon spin 
polarization P z (t), which can be detected by counting 
emitted decay positrons forward (f) and backward (b) 
of the initial muon spin direction; this is possible due 
to the asymmetric nature of the muon decay [13, flij |. 
which takes place in a mean time of 2.2 /is. In our 
experiments positrons are detected by using scintilla- 
tion counters placed in front of and behind the sam- 
ple. We record the number of positrons detected by for- 
ward (Nf) and backward (iVb) counters as a function of 
time and calculate the asymmetry function, G z (t) , using 
G„(t) = [N f (t) - a cxp N h (t)]/[N { {t) + a exp N b (t)], where 
a oxp is an experimental calibration constant and differs 
from unity due to non-uniform detector efficiency. The 
quantity G z (t) is then proportional to P z (t). 

Raw data for EuO at two temperatures, one slightly 
higher and one much lower than Tq, are shown in 
Fig. Ufa). These demonstrate that spin relaxation above 
Tq changes to a damped coherent oscillation below Tq . If 
the experiment is repeated with a sample of EuO in which 
0.6% of the Eu ions are replaced with Gd, an almost iden- 
tical result is obtained [Fig. [ljb)] although the damping 
rate of the oscillations in the ordered state is noticeably 
larger (by about a factor of two). This level of Gd dop- 
ing is known to introduce a static magnetic inhomogene- 
ity [lj|. The frequency of the oscillations in both cases 
increases rapidly on cooling below Tc, see Fig. [TJc,d,e), 
approaching w 30 MHz at zero temperature. The Fast 
Fourier transform data in Fig. [He) demonstrate that only 
a single precession frequency is observed. The tempera- 
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FIG. 1: (Color online.) (a) Raw ^SR data for EuO. Above T c 
(at 70 K), -P 2 (t) relaxes. An oscillating signal develops below 
Tc (data shown for 1.5 K). (b) The same for Euo.994Gdo.006O. 
(c) Plot showing FFT of the muon data for EuO as a function 
of temperature. The single precession frequency which follows 
the order parameter is clearly visible, (d) The extracted pre- 
cession frequency for EuO as a function of temperature, (e) 
The same experiment repeated for Euo.9g4Gdo.006O. The line 
through the data in (d) and (e) is identical for the two samples 
and is a best fit of the data to the phenomelogical function 
v(T) = i/(0)(l - (T/Tc) a ) 13 , producing a « 1.5, j3 w 0.4). A 
more reliable extraction of the critical exponent ft, focussing 
only on the critical regime, is given later in the paper. 



ture dependence of the precession frequency in the Gd- 
doped sample is almost identical to that in the pure sam- 
ple, demonstrating that the order parameter is following 
the intrinsic magnetism in the EuO host and is relatively 
insensitive to low levels of doping. In fact, at the 0.6% 
level of Gd-doping, 93% of the Eu ions have their full 
complement (twelve) of nearest-neighbour Eu ions, with 
most of the remaining Eu ions having only one of those 
nearest neighbours replaced by Gd. For both samples, 
the relaxation rate of the oscillatory component rises as 
T approaches Tq from below and the low-temperature re- 
laxation rate of the Gd-doped sample is larger than that 
of the pure sample 

The muon spin precesses around a local magnetic field, 
(with a frequency v — ( i y li /2'ir)\B tl \, where 7^/271 = 
135.5 MHz T^ 1 ). This local field (B^ w 0.22 T at T = 0) 
is a sum of various terms, including the Lorentz field £>l, 
the hyperfine field £?hf, the demagnetizing field -Bdemag 
and the dipolar field Bdi p . The latter quantity is a func- 



tion of the muon-site r M and can be written as 



B£ p (iv) = £A ap (*v)™ 



a/3. 



(1) 



a sum over the magnetic ions in the crystal; the magnetic 
moment of the ith ion is m^. In Eq. ([1]), D"^(r M ) is the 

dipolar tensor given by I)f (rj = j^jt (^%t^- - , 
where Ri = (Rf,Rf,R*) = r M — r^. The behaviour of this 
tensor is dominated by the arrangement of the nearest- 
neighbour magnetic ions and leads to a non-zero local 
magnetic field for almost all possible muon sites [20j]. On 
electrostatic grounds, the likely muon site in EuO is at 
the ±11 position [see Fig. Ufa)], equidistant from four 
Eu cations and four oxygen anions. Because of the mag- 
the easy axis for the Eu moments 



21 



netic amsotropy 
is along the (111) set of directions and for this moment 
alignment the \ \\ position is a point at which the dipo- 
lar magnetic field actually vanishes. The value of i?di P 
has been calculated for the case in which the muon is 
displaced from the \ \\ site and its position is allowed to 
vary along the [111] direction, see Fig.[2jb). The dipolar 



field vanishes at both \\\ 



[muon site) and \\\ ( ox y _ 



gen site) and increases sharply as the site moves away 
from these special positions of high symmetry. The two 
curves show the cases in which the muon displacement is 
the same or a different choice of (111) direction. In the 
former case, the dipolar field at the muon site is parallel 
to the moment direction; in the latter case (for which 
there are three possibilities), it lies along one of the crys- 
tallographic axes and its amplitude is reduced by a factor 
of -\/3- Thus if the muon site is displaced from the \\\ 
position towards a particular oxygen anion, then there 
would be a contribution to the dipole field resulting in 
two precession frequencies with a ratio of v3j and with 
amplitudes in a 3:1 ratio. Since only a single frequency is 
observed, we conclude that the muon site is indeed at the 
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FIG. 2: (Color online.) (a) Muon site in EuO. The crystal 
structure for one-eighth of the unit cell is shown (the side of 
the cube is a/2), (b) Calculated dipole field along the (111) 
axes in EuO. The muon is located at £££ and the Eu moments 
are all aligned along [111] (solid black line) or one of [111], 
[111] or [111] (dashed red line). 



\\\ position in which the dipolar field is zero, so that 
the observed local field (B^ = 0.22 T at T = 0) is due 
to a sum of the Lorentz field (£?l = = 0.80 T at 

T = 0), -Bdemag and £?hf- Since the sample is polycrys- 
talline and multidomain, we neglect -Bdemag and deduce 
that the hyperfine field -Bhf < (antiparallel to the mag- 
netization), as found for EuS [22j, and takes the value 
-Bhf = -B L ± B^, and so either -0.58 or -1.02 T. For 
both samples, the amplitude of the oscillatory compo- 
nent is reduced from the full value at low temperature 
[see Fig.[lja,b)], but recovers on warming towards Tc, so 
a fraction of muons may implant in some additional state 
which depolarizes the muon very rapidly. 

We note that a recent experiment 23] on SmS has 
shown evidence for the formation of a bound magnetic 
polaron consisting of an electron around the implanted 
muon, in which the electron localization is stabilised by 
exchange energy. This occurs in the paramagnetic state 
in which a ferromagnetic droplet is localized in the para- 
magnetic host. A similar effect has been noted in EuS 
241 ] although the larger magnetization ensures it occurs 
at temperatures 3> Tq ^25j and the same will be true in 
EuO in which the magnetization is even larger. There- 
fore such a muon-related polaron is not relevant for EuO 
in the studied temperature regime. 

The temperature dependence of the precession fre- 
quency for both EuO and Euo.9g4Gdo.006O was followed 
near Tc and the results are plotted in Fig. [3] The fitted 
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FIG. 3: Precession frequency extracted from /iSR data as 
a function of temperature close to 7c for (a) EuO and (b) 
Euo . 994 Gdo .006 O . 



values of Tc and the critical exponent /3 are similar in 
each case, though the value of ft is quite sensitive to the 
precise value taken for T c . Due to the difficulty in sta- 
bilising the temperature better than ss 10 mK, we do not 
believe that the difference between the two values of j3 
is significant. They are both close to 0.36-0.37 obtained 
using neutron scattering |26| and 0.38 obtained from a 
second order e expansion for the Heisenberg ferromagnet 
with dipolar interactions [271 ] . 

Above Tc, we observe simple exponential relaxation 
[Fig. [T{a,b)] with a relaxation rate A. For zero-field re- 
laxation of muons initially polarized parallel to z, A can 
be written in terms of field-field correlation functions us- 
ing A = 4-1-00^ ((B x (0)B x (t)) + (B y (0)B y (t))) M- 
When each Eu spin component fluctuates, it produces a 
field fluctuation via the resulting modulation of the dipo- 
lar and hyperfine couplings. Our measurements of the 
zero-field relaxation rate for the EuO sample are plotted 
in Fig. 2] (a less complete set of data for the Gd-doped 
sample is also shown). There is a small rise in A as the 
temperature is lowered towards Tc but apart from this A 
remains just below sa 2 MHz for both samples across the 
entire range studied. 

These results can be compared with calculations on a 
localized Heisenberg ferromagnet which have been per- 
formed with EuO in mind [la, 16] (Fig. 4). The theory 
of Lovesey and Engdahl [l5j includes only the dipolar 
coupling and has been evaluated for temperatures above 



1.3Tc, assuming a muon site of \\\- Though underes- 
timating the observed experimental values, this theory 
does remarkably well in providing a good estimate of the 
rough size of A, the discrepancy perhaps being due to 
neglecting the hyperfine contribution. It is known that 
critical fluctuations enhance the role of the hyperfine cou- 
pling over the dipole coupling [14[ because -Bdip = at 
the muon site in the ordered state and the peak in the 
susceptibility is at k = 0. Nevertheless, when the dipo- 
lar calculation is extended into the critical regime (just 
above, and very close to, Tc) it predicts a divergence in 
A which is not observed. An earlier mode-coupling ap- 
proach [l6| also predicts a very sharp increase in A on 
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FIG. 4: (Color online.) Relaxation rate as a function of tem- 
perature in the paramagnetic regime for pure EuO (solid cir- 
cles) and (b) Euo.994Gdo.006O (dots). The predicted relax- 
ation rates for T > Tc according to Ref. LL5] (numerical cal- 
culation: thick solid line, purple; critical regime using experi- 
mental values of correlation length: thick dashed line, purple) 
and Ref. [l6| (dotted line, blue) are also shown. 



cooling to Tc from about 0.3 K above it; this is also not 
observed in our data. Magnetic polaron formation has 
been detected using Raman scattering [28| in a narrow 
range (~ 20 K) above To It may be that the forma- 
tion of magnetic polarons modifies the relaxation in this 
regime from that which would be expected from theory, 
perhaps by providing an additional relaxation channel 
for the muon which masks the critical slowing down pre- 
dicted by the theory and hence the absence of the diver- 
gence in A. We note that a similar absence of a divergence 
in A is observed in EuBg 29] in which magnetic polarons 
have been found [jjjj . 

In conclusion, we have identified the muon site in EuO 
and estimated the hyperfine field. Our results confirm 
long-range order which is relatively insensitive to low 
doping of Gd. The measured A in the paramagnetic state 
agrees quite well with the theory of Lovesey and Engdahl, 
but the available theories fail in the critical regime, pos- 
sibly due to magnetic polaron formation. 
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